Doping dependence of the Neel temperature in Mott-Hubbard antiferromagnets: 

Effect of vortices 
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The rapid destruction of long-range antiferromagnetic order upon doping of Mott-Hubbard antifer- 
romagnetic insulators is studied within a generalized Berezinskii-Kosterlitz-Thouless renormalization 
group theory in accordance with recent calculations suggesting that holes dress with vortices. We 
calculate the doping-dependent Neel temperature in good agreement with experiments for high-Tc 
cuprates. Interestingly, the critical doping where long-range order vanishes at zero temperature is 
predicted to be Xc ~ 0.02, independently of any energy scales of the system. 
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The study of lightly doped Mott-Hubbard (MH) anti- 
ferromagnetic insulators is of great current interest, since 
the insulating parent compounds of cuprate high-Tc su- 
perconductors are of this type. The various parts of the 
phase diagram of these compounds are believed to be in- 
timately related. Therefore it is important to understand 
the properties of the antiferromagnetic phase, in particu- 
lar the rapid destruction of magnetic order upon doping 
and the anomalously small critical doping of Xc ~ 0.02 
holes per copper ion |i|. In the present letter we de- 
rive the Neel temperature T/v as a function of doping us- 
ing a generalized Berezinskii-Kosterlitz-Thouless (BKT) 
renormalization group theory MM for the vortices in the 
antiferromagnetic state. There are two types of vor- 
tices, thermally created electrically neutral ones and elec- 
trically charged ones, which are centered at the holes. 
Both types are nucleated separately (as vortex-antivortex 
pairs), but additively screen the vortex interaction, with 
a common unbinding temperature Tj^{x). This temper- 
ature is indeed strongly reduced upon doping and van- 
ishes at Xc ~ 0.02 independently of the energy scales of 
the system. Our approach is independent of any partic- 
ular microscopic model and can thus serve as a guide for 
electronic theories. 

Our physical picture is the following: The holes in- 
troduced by doping are mainly located at the planar 
oxygen sites, where they frustrate the antiferromag- 
netic exchange interaction between copper spins due to 
their tendency to form copper-oxygen spin singlets [Q,^ . 
This may lead to ferromagnetic coupling between the 
two spins l^]. Since the system is approximately two- 
dimensional, the staggered magnetization can form a vor- 
tex as sketched in Fig. |l| to evenly distribute the frustra- 
tion induced by the hole. On the other hand, neutral 
vortices without a hole in their core can be thermally 
created. Due to the easy-plane Dzyaloshinskii-Moriya 
anisotropy m, the staggered magnetization can be de- 
scribed by a two-component order parameter at low en- 
ergies, leading to a logarithmic size dependence of the 
single- vortex energy P,P[. This implies a logarithmic 
vortex interaction, making BKT scaling ideas applica- 



ble. To describe the interplay of charged and neutral 
vortices, which determines the Neel temperature Tjv, we 
have to extend the BKT theory to a system with two 
kinds of topological defects. One important point is that 
the density of charged vortices is given by the doping x. 
The other is that we have to describe the screening of the 
vortex interaction due to both types. 
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FIG. 1. Schematic representation of a charged vortex in- 
duced by a single hole at an oxygen site (black circle). The 
hole frustrates the interaction between copper spins (arrows). 
The circles denote oxygen atoms. A neutral vortex would be 
centered on the square between four oxygen atoms and not 
have a hole in its core. 

Electronic theory supports our picture: Within unre- 
stricted Hartree-Fock theory Verges et al. [g[ found sev- 
eral competing low-energy configurations for a lightly 
hole-doped Hubbard model, including spin polarons, do- 
main walls, and holes dressed with vortices and antivor- 
tices. In a spin polaron the staggered local moments in 
the vicinity of the hole are reduced but still coUinear, 
while in a vortex (antivortex) they rotate through 27r 
(— 27r). Seibold pU using a slave-boson approach and 
Berciu and John loTwithin a self-consistent Hartree-Fock 
theory found that an even number of holes dress with vor- 
tices and antivortices (or merons P| ) in the ground state 
for appropriate parameters. Since the energy of a single 
vortex diverges (logarithmically P,H,Hl) with system size, 
whereas that of a vortex-antivortex pair remains finite. 



only pairs are created in infinite systems. 

An advantage of BKT-type tfieories is that they in- 
clude fluctuations on all length scales, in particular on 
large ones, which are crucial close to the phase transition. 
Previous studies that only included fluctuating local mo- 
ments without correlations between them overestimated 
the critical doping Xc uM ■ Similarly, the decrease of T^ 
with X obtained from the fluctuation exchange approx- 
imation is also too slow |12]. By including correlations 
between neighboring spins into a slave-boson theory for 
the three-band Hubbard model, Schmalian et al. ob- 
tained a critical doping of Xc ^ 0.025 in better agreement 
with experiment. However, this approach takes only fluc- 
tuations on the length scale of the lattice constant into 
account, but neglects fluctuations on larger length scales. 

In BKT theory the interaction of vortices is screened 
by the polarization of vortex pairs lying between them 
[g,||. As the temperature is increased more pairs are 
thermally created leading to increased screening. At the 
Ncel temperature the screening becomes strong enough 
for the largest pairs to break up. The resulting free vor- 
tices destroy the magnetic order. The situation here is 
more complex: First, we have to take account of the 
screening due to both charged and neutral vortices, and 
second, the density of charged vortices is fixed by the 
doping level x. 

Only neutral pairs are thermally created, whereas 
charged vortices enter the system only upon doping. In 
principle one could imagine a single hole doped into 
the system to form a pair consisting of a charged vor- 
tex and a neutral antivortex or vice versa, but micro- 
scopic calculations |P,0,p| do not find this configuration 
at T = 0. Rather, two holes are needed to produce a 
vortex-antivortex pair. For simplicity we assume this to 
hold also at finite temperatures. However, even if mixed 
pairs are not created upon doping, they are formed when 
vortex pairs exchange partners. 

The density of neutral vortices is controlled by their 
chemical potential /inou or equivalently the vortex core 
energy E'coro = — /^neu, which depends on details of 
the copper-oxygen and copper-copper interactions and 
is treated here as a parameter. The energy of a vortex- 
antivortex pair is 2i?corc + V with the interaction M] 
V{r) — q^ ln(r/ro), where q^ = 2ttJS^ is the strength 
of the interaction, J is the exchange interaction between 
nearest neighbors, S = 1/2 is the spin, and ro is the 
small-distance cutoff of BKT theory, ro can be inter- 
preted as the smallest possible vortex-antivortex separa- 
tion PJq]. Two charged vortices additionally experience 
a Coulomb interaction, which, however, is irrelevant in 
the renormalization-group sense, since it falls off faster 
than V{r) ^ Inr. 

The probability of creating a neutral or charged vor- 
tex in an area Tq is given by its fugacity y^cu and ych, 
respectively. Since we have assumed that vortices are 
only created in neutral or charged pairs, we consider the 



pair fugacities y^^^ and y'^^. For the smallest possible 
neutral pairs of size rp. 



y2 (r^) ^ C^ e^''^"™ 



(1) 



where Cncu is a constant of order unity [^ and (3 is 
the inverse temperature. The constraint on the density 
of charged vortices is implemented by choosing ychl'^o) 
in such a way that their total density equals the hole 
density, see below. The vortex interaction is screened 
by the polarization of smaller vortex pairs, V{r) — 
q^ / e{r) ln(r/ro). The screening is described by the spin- 
wave stiffness K{r) — f3q'^/2TTe{r). In the renormaliza- 
tion group, small pairs of sizes between r and r + dr are 
integrated out and their effect is approximately incor- 
porated into renormalized quantities K{r), j/ncu(^)' ^^^^ 
y%Jr). Starting from r — rp, this operation is repeated 
for larger and larger pairs leading to the recursion rela- 
tions i,|ll 



dy^/dl 
dK/dl 



2i2-7:K)yl 



(2) 
(3) 



The initial conditions are Eq. (m and K{rQ) — (3q^ /2'k. 
Equation (0) determines the fugacities of neutral (77 — 
neu) and charged (77 = ch) pairs of size r = roe'. Two 
separate equations for y^^^ and y^^i ^^^ present, since we 
assume that vortices are created either as neutral pairs or 
as charged pairs with two holes. Both types feel the same 
screened interaction V at large distances so that the same 
stiffness K appears. Differences at smaller separation 
are incorporated into the core energies. Equation (|^) 
describes the additional screening due to pairs of size 
roe'. Their total density is proportional to y^^^ + y^^i- 

If the stiffness K vanishes for Z ^ 00, the interaction 
is fully screened for large pairs (e — > 00), which thus be- 
come unbound, destroying the magnetic order. Since the 
interaction on large length scales is the same for neu- 
tral and charged vortices, this unbinding happens at a 
single transition for both types. While solving Eqs. (||) 
and (H) we have to simultaneously satisfy the constraint 
on the density rich of charged vortex pairs. As shown in 
Ref. |13], this density can be expressed in terms of the 
fugacity. 
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dle-''yt^{l). (4) 



The pair density has to equal half the density of holes, 
"■ch = x/2a?' , where a is the lattice constant. In practice, 
the recursion relations are integrated numerically to find 
y%Jl), from which we calculate rich. The initial value 
j/^j^(0) is varied until the contraint is satisfied. 

The resulting phase diagram is shown in Fig. 0. We 
used Cnou = 1, J^ = 1800 K, and ro — 2a, and varied 
the core energy -Ecoro- The phase below the Neel temper- 
ature T/v shows quasi-long-range antiferromagnetic or- 
der, which is made long range by the weak interlayer 



exchange. The phase for T > T/v or x > Xc is char- 
acterized by free vortices, which destroy the long-range 
order, but leave short-range order on the length scale of 
the separation between free vortices intact. Short range 
correlations have indeed been observed in cuprates up to 
much larger dopings. For T — > we find the critical pair 
density, 



^,=h« 0.04273 V 



(5) 



The numerical factor is universal: Since neutral vortices 
do not exist for T — > 0, it cannot depend on -Ecorc and 
Cncu- The remaining energy scale q^ does not enter the 
result, since it is multiplied by the diverging /3 — l/ksT. 
While n^j^ and thus the critical doping Xc = 2n'^^a^ are 
independent of energy scales, they do depend on the non- 
universal minimal separation rg of two vortices which is 
of the order of twice the core radius or twice the corre- 
lation length. Slave-boson and Hartree-Fock calculations 
[jlO|,P[ show that the core radius is not significantly larger 
than a lattice spacing a. In Fig. |2| we have taken the 
core radius to be a so that tq = 2a, which results in 
Xc ~ 0.021 in very good agreement with experiments on 
La2-xSr^Cu04 ||,|l6|-|l. For Yi^^Ca^BaaCuaOg exper- 
iments find Zc/2 ~ 0.03 [|8[, of the same order of mag- 
nitude as our value of 0.021 (the number of holes per 
copper atom is z/2 in this double-layer compound). For 
YBa2Cu307_i- a critical hole concentration of 0.021 cor- 
responds to Sc ^ 0.68, following Tallon et al. [M, in good 
agreement with experiments pQ]. 



500 






400 



300 : 



g_ 200 

E 

100 







)c//4 






\ 




£c„re = ?/2 


\ 

\ 




£c„re = ?/4 


-, \ 






% \ 










{q^ = 2njf) 


- \ 




paramagnet 




\ T 


; 


antiferromagnet 


\ 






X , , 







0.01 



0.03 



0.04 



0.02 
doping X 

FIG. 2. Phase diagram for the antiferromagnetic phase of 
the cuprates. The Neel temperature Tjv is calculated using 
four values of the core energy -Ecore of thermally created vor- 
tices. Note that the curves for larger E^oic are very close to 
the one for Ecoie = 2q^. The symbol x at zero doping denotes 
the maximal possible transition temperature. 

In the high-temperature region of the phase diagram 
the overall temperature scale is set by the maximal pos- 
sible Neel temperature T^'^'^ = q^/A = nJS^/2. For 
S = 1/2 this gives T^'''' « 0.393 J, compared to the 



mean-field result T™^ = 0.5 J for a Heisenberg anti- 
ferromagnet on a cubic lattice with interlayer exchange 
J± <C J. The reduction of T/v is due to fluctuations, 
which are strong for quasi-two-dimensional systems. The 
actual value of T/v (a; = 0) and the shape of the curve 
T/v (a:) at small doping are determined by the core en- 
ergy Ecoic- Note, we obtain the correct temperature scale 
under the reasonable assumption that i^core is not very 
much smaller than the interaction strength g^. Experi- 
mentally, T/v is found to depend only weakly on doping 
for small x pJ\ ^ , which requires a small core energy. Then 
many neutral vortices are created at a given temperature 
so that charged vortices only become relevant at higher 
doping. Conversely, for large core energy only a few ther- 
mal vortices are present even at T^'^^. For i?core ^ 2q'^ 
the curve Tn{x) in Fig. g does not change appreciably 
with i^corc so that the charged vortices would determine 
the magnetic properties even at very small doping. 
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3. Comparison of the calculated Neel temperature for 
EcoTc/q^ ~ 0.13 (solid line) and results for Tjv from NQR 
and /iSR experiments for La2-a;Sra:Cu04 (triangles) |l7|. The 
dashed line schematically shows the freezing temperature Tf 
[hsl, below which the holes become immobile. 

Quantitative agreement with /iSR and NQR experi- 
ments on La2-a:Sra;Cu04 by Borsa et al. p7| can be ob- 
tained by appropriate choices of the exchange J, the 
core energy Ecorc, and the core radius tq, see Fig. [|. 
For this plot, J = 2410 K, Score = 0.13039^ = 493 K, 
and ro = 2.052a. Typical experimental values are J ^ 
1400 K pl|. This discrepancy may be due to the sim- 
plified description of the anisotropic antiferromagnet by 
two-component spins or to the neglect of the interlayer 
exchange J± and the doping dependence of J. 

We now briefly comment on electron-doped cuprates. 
While we reproduce the correct critical doping for hole- 
doped compounds, our approach would not yield the 
much larger critical doping Xc,e ~ 0.14 in electron-doped 
cuprates ||l|. The reason is that the additional electrons 
mainly fill up the copper 3d orbitals and destroy the mag- 



netic moments at the copper sites. Thus, the main effect 
of electron doping is to dilute the antiferromagnet. There 
is no spin-singlet formation involved and hence no ten- 
dency towards vortex formation. Then differently as in 
the case of hole doping T/v decreases due to spin dilution. 

The validity of our approach is questionable if the 
charged vortices become immobile. There is experimen- 
tal evidence that this happens below a crossover tem- 
perature Tf, which increases with doping and reaches 
about 16 K in La2-i:Srj^Cu04 for x w Xc and falls off 
again for larger x [p2|,^. Tj is sketched in Fig. 0. Be- 
low this temperature the holes (charged vortices) form 
a glass and their dynamics strongly slows down. The 
experimental Neel temperature in this region should de- 
pend strongly on the time scale of the experiment. On 
the other hand, the true phase transition is governed by 
the behavior in the limit of infinite time so that the for- 
mation of a glass below Tf should affect it only weakly. 
The doping dependence of Tf can be qualitatively un- 
derstood in our picture: In the magnetically disordered 
phase the logarithmic part of the interaction of charged 
vortices is screened on the length scale of the correlation 
length ^,g,|l5[, which is still large close to T/v. Thus the 
interaction of charged vortices (holes) changes smoothly 
at the transition and decreases for larger doping, leading 
to similar behavior of the freezing temperature Tf . 

There are theoretical indications that the holes may 
form one-dimensional stripes at low temperatures [GSJ. 
Modeling stripes by a phenomenological anisotropic Hei- 
senberg model, Castro Neto and Hone [|^ calculated 
the Neel temperature Tn(x) within a renormalization- 
group scheme and found good agreement with experi- 
ment. However, in this theory the critical concentration 
Xc is basically a free parameter. Note, stripes formed by 
charged vortices consist of alternating vortices and an- 
tivortices, in order to lower the interaction energy. These 
stripes are automatically anti-phase domain walls P5| , 
which are observed experimentally. 

In conclusion, by starting from the assumption that 
holes doped into the Mott-Hubbard antiferromagnet 
dress with vortices and using independently obtained val- 
ues for the exchange interaction and the antifcrromag- 
netic correlation length in the ordered phase, we obtain a 
doping-temperature phase diagram for the antiferromag- 
netic phase in qualitative agreement with experiment. In 
particular, the predictions for the critical doping at zero 
temperature and the Neel temperature at zero doping 
are of the observed order of magnitude. Our approach 
uses a generalized BKT theory, which does not depend 
on any particular microscopic model. With an appropri- 
ate choice of the core energy of thermally created vor- 
tices we can obtain quantitative agreement with experi- 
ment. The core energy controls the shape of the phase 
boundary at small doping, but does not affect the re- 
gion of higher doping, where the Neel temperature ap- 
proaches zero at a critical hole density that is univer- 



sal in natural units. Our results show that stripes are 
not required to understand the data. The success of 
this theory based on vortex fluctuations emphasizes the 
importance of two-dimensionality for understanding the 
cuprates. This should also hold in the more strongly 
doped superconducting region. It would be desirable to 
include the spin rearrangement around holes, which is 
induced by frustration due to singlet formation, into an 
electronic theory of underdoped cuprates. For the an- 
tiferromagnetic phase such a theory should yield results 
similar to the ones shown here. 
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